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A method is described for obtaining the 
rate constant of the photodegradation 
process of fluorophores illuminated 
by a focused laser beam. The expUcit 
kinetic equations, describing the 
population dynamics of excited singlet 
and triplet states, are averaged over the 
illuminated volume to describe the 
resulting fluorescence signal. The 
illumination is modulated at frequencies 
from 1 Hz to 100 Hz. Synchronous 
detection of the resulting fluorescence 
yields in-phase and quadrature compo- 
nents. The measurement of the ratio of 
quadrature to in-phase components at 
several power levels yields information 
on the photodegradation rate. Specifically 
it is shown that the data can be interpreted 
in a manner which yields the value of 
the photodegradation rate independently 
of other parameters entering the model. 
Experiments are performed with erythro- 
sine B which has a large intersystem 
crossing rate to the triplet state. 
Measurements in solutions with different 
viscosities show that the photodegradation 
rate depends on the viscosity. This is 
interpreted as evidence for an inter- 
molecular interaction mechanism. We 



explore the uncertainty of the estimated 
photodegradation constant taking into 
account the uncertainties of the measure- 
ments used in the synchronous detection 
technique. 
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1. Introduction 

This work focuses on the photodegradation of fluores- 
cent fluorophores in solution illuminated by a focused 
laser beam. The technique described in this paper may 
provide a quantitative method for studying laser-induced 
effects. There are many instances of a laser beam inter- 
acting with biological systems [1]. Examples range from 
laser surgery to the simple probing of labeled antibodies. 
In most cases the interpretation of the laser-based meas- 
urements is rather difficult. Fluorescent labels are wide- 
ly used to probe cell structure and cellular events, and 
photodegradation is often a limiting factor There are 



many difficulties in interpreting the reduction of fluores- 
cence signal due to photodegradation. Heating, convec- 
tion, and scattering are some of the confounding factors. 
We describe a method to obtain quantitative information 
about the photodegradation rate constant of fluo- 
rophores. Another example of an application that may 
benefit from this methodology is photodynamic therapy 
[2]. Our study extends the work of a previous paper [3] 
which presented a derivation of the kinetic equations 
describing photochemical reactions in a focused laser 
beam. This paper shows how harmonic solutions of the 
kinetic equation can be used to analyze experimental 
measurements and extract photodegradation rates. 
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2. Formulation of the Model 

Figure 1 shows a schematic of the physical situation 
under consideration in this work. A more detailed 
description of the apparatus is given in reference [4]. A 
laser beam, incident along the z axis, is focused to a waist 
of about 20 |am. A solution of fluorophores flows along 
the X axis (vertical direction in Fig. 1 ) and passes through 
the focused laser beam. A detector (not shown in Fig. 1) 
located along the y axis, which is perpendicular to the 
plane of the figure, detects fluorescence from the illumi- 
nated fluorophores within the acceptance aperture of the 
detector The acceptance aperture should be no larger 
than twice the Raleigh range defined as the distance over 
which the beam diameter does not exceed 1.414 of the 
value of the minimum laser beam diameter (called the 
beam waist). The detector sums the photons coming 



from the fluorophores in the observation region. The 
kinetic equations that describe the temporal evolution of 
the fluorophore concentration are given for a three state 
model by 

dN,{r,t) ^^ dN,{r,t) _ 
dt dx 

-k^(r,t)N^(r,t) + k^N,(r,t) +k,N,(r,i) , 



dt dx 

+K{r,t)N,{r,t)-k^N,{r, t)-k,,^N,{r, t) , 

dt dx 

+k,^^N,ir,t)-k,N,ir,t)-k,N^ir,t). 



(1) 
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Fig. 1. A schematic of the experimental arrangement described by the analysis. A laser 
beam propagates in the z direction (to the right) and is focused to a narrow waist. The flu- 
orophore solution flows in the x direction (vertical in the figure). The detector is located 
along the>' axis (perpendicular to the plane of the figure). The fluorescence is measured 
from a portion of the illuminated region as defined by the detector aperture. 
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The symbol r represents the x,y coordinates, and v is 
the flow velocity. The various rates in Eq. (1) are defined 
in Fig. 2 which shows a schematic of the three state 
model under consideration. The radiative absorption rate 
is given by k^, while the total relaxation rate to the 
ground state is given by k,. The symbol A;,^^ represents the 
transition rate between the singlet and triplet excited 
states. The radiative and non radiative relaxation of the 
triplet state are represented by kj- and kj respectively. The 
boundary conditions are given by 

7V„ i-L, y, z, t) = N\ N, i-L, y, z, t) = 0, 

N^{-L,y,z,t)=0 for t>0 
N, (x, y, z, t) = N\ N^ (x, y, z, t) = 0, 

TVj {x, y, z,t)=0 for x<-L 

N,ix,y,z,0)=0, N,i,x,y,z,0)=0 
for x>-L 

The first two lines of Eq. (la) state that prior to enter- 
ing the laser beam located at x = 0, the fluorophores are 
in the ground state with a concentration given by N". 
The third line of Eq. (la) states that prior to turning on 
the laser at / = 0, the excited state populations are zero. L 
is a distance beyond which the laser intensity vanish- 
es. The inclusion of the triplet state provides a more real- 
istic model of the photodegradation process than the two 



(la) 



state model. The absorption rate, k^, shown in Fig. 2 is 
the only rate constant which varies with time and posi- 
tion. The solution of the extended model will follow that 
given previously for the two state model [3]. Briefly, 
Eq. (1) is multiplied by a constant s= 10"* and it is 
noted that the factors k,s and kj,^s are of the order 1, 
while kjr,t)s is of the order 1 over most of the illumi- 
nated volume; kj€ is of the order 10"^ and k^s is of the 
order 10"^. Therefore it is convenient to look for solu- 
tions to Eq. (1) which are correct to different orders in s. 
Specifically, we write the solutions to Eq. (1) as 

Noix,y,t,s) = N^ix, y, t) +eN^i{x, y, i) + 0{e^) 

N, (x, y, t, e) = N, (x, y,t)+sN,^ (x, y,t) + 0(s^) (2) 

N, (x, y, t, s) = N,ix,y,t)+sN^iix,y,i)+0{s^). 

Inserting Eq. (2) into Eq. (1) and collecting different 
orders of s, we obtain for the zero order the relations 
shown as 



N,ir,t)=^N,ir,t) 

N,ir,t) = ''f/-"l-^ N,ir,t). 



(3) 



N, 




N 







Fig. 2. The molecular energy levels involved in the photochemical process. A^q, N^, and A'3 are the populations of the ground, the lowest excited 
singlet, and the lowest excited triplet states of the molecule. The arrows indicate the possible transitions between the states and the letters next to 
the arrows represent the rates of the transitions. The thick line extending the excited singlet state to the right represents the manifold of states 
involved in the equilibration of the initial excited singlet state. The rate k^ represents all processes that result in a non fluorescent fluorophore. 
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Next, the total number of intact fluorophores is 
defined by 



N{r, t) ^ N, (r, t) + N, (r, f) + N,{r,t). 



(4) 



Thus to the first order in £, the kinetic equation for 
N{r,i) is given by 



dN{r,t) dN(r,i) , - , , 



3/ 



dx 



fluorophores integrated over the entire observation 
region, {N{t)), is reproduced below [3]. 



y,,) 



-riajP^+Ajt)) 
l + ba(Pg+A(t)) 



{Nit)) 



V (l + baXP,+A(t))) 
1.274 -w (l + ba(P^+A(t))) 



(9) 

i{Nit))-N') 



where 



Eq. (3) and Eq. (4) are used to rewrite N^(r,t)m terms 
of N(r, i) and the resulting expression for N^ (r, i) 
is inserted into Eq. (5). The result is given by E q. (6) 
where in addition we assumed that kj. <c k^^, and terms 
of the order k^ have been dropped. 



dN{r,t) dN{r,t) 



dt 



dx 



Kir,t)K 



(6) 



' K{K+KX\ + ^4r^) 



N{r,t). 



kj. {k^ + ki^^ ) 



Using the same notation as in reference [3] define 



kSr,t) = (yP^P(r,t) 



k^(l + ^-) 



(J) 



ri = k,b 



Where P{r,t) is the laser power, a (cm ) is the molec- 
ular absorption cross section, P^. = Xn/hc where X is the 
wavelength, n is the index of refraction, h is Planck 
constant, and c is the speed of light in vacuum. 
Inserting the above definitions into Eq. (6), we obtain 
the final form of the kinetic equation. 



dN(r,t) dN(r,t) 



dt 



dx 



'\ + bP{r,t) 



N{r,t) (8) 



A(/) = /J cos(ft)0 and ( A^( ) = fj /( X y) N( x, y, t) dxdy 

R 

Here Pi is the amplitude of modulated laser power, Pq 
is the constant laser power, /(x,>') is the spatial profile 
of the laser beam, and co=2 nf, the frequency of 
modulation. The other terms in Eq. (9) are V , the veloc- 
ity of flow, w, the width of the laser beam, and b, and 
r\ are defined in Eq. (7). The quantities a and a 'depend 
on the beam profile function/(x,>') which is assumed to 
be a Gaussian function with area normalized to 1. 
Explicitly 



a = 



J/(-^> j)/(-^> y)N{x, y)dxdy 
jf(x,y)N(x,y)dxdy 



(10) 



The value a' a of is calculated using Eq. (10) with 
N{x,y) replaced by the derivative dN{x,yydx. It was 
shown in a previous paper [3] that to a very good approx- 
imation the values of a and of are time independent 
when r\ is small, and that they can be calculated using the 
time independent, steady state solutions of the complete 
set of kinetic equations. Consequently we use the time 
independent solution, N{x,y), in Eq. (10). The spatial 
profile of the laser beam is described by the function [5] 



.. X 1 . x'-+y\ 
f{x,y) = ^exp( T— ) 

nw w 



P^{x,y) = PJ{x,y) 
AP{x, y, t) = Pif{x, y) cos{cot) . 



(11) 



Equation (8) is identical to Eq. (14) in Ref [3]. The 
averaging of Eq. (8) over the observation volume 
follows the same procedure as shown previously [3]. 
The resulting kinetic equation for the number of intact 



The fluorescence signal originates from the radiative 
relaxation of molecules excited by the laser beam. The 
fluorescence signal from any location in the laser beam 
depends on the laser power and the concentration of 
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fluorophores at that location. An explicit form of the 
fluorescence signal integrated over the entire observa- 
tion region can be derived by an averaging argument. 
The result is given by Eq. (12) and is accurate for small 
laser power levels. 

F{i) = B%^,b[p{t){N{t))-ba{P{t)f{N{t))\. (12) 

Here S'is an instrument constant, k^^^ is the radiative 
decay rate of the molecular excited state and P{t) = 
Po+ P\ cos {(Ot) =Po + A(i). The fluorescent signal can 
therefore be obtained (with a small error) by solving a 
single ordinary differential equation, Eq. (9), for {N{t)), 
and then using the solution in Eq. (12). Together these 
equations represent an experimentally accessible math- 
ematical model of the fluorescent signal due to excita- 
tion by a focused laser beam. The model can be used to 
analyze measurements. The problem is that the two 
equations contain many parameters, r\, b, a, a', v, and 
w, and it is not clear how to extract the value of the 
relevant physical quantity, k^. 



3. Proposed Measurement Procedure 

3a. Derivation of the Harmonic Measurement Model 

Consider the case of very small laser power We make 
the approximation, baP(i) « 1, ba'P(t) « 1, and write 
Eq. (9) as 



im^ 



riaiP,+P,e>'") + 
{Nit)) = 



1.274 -w 



1.274-w. 



(13) 



Where j is the imaginary number, and we have 
written the modulated laser power as a complex quan- 
tity, P(t) = Po + Pye'"". By writing the time dependent 
power modulation as P^e''"' (Py real) and assuming a 
complex solution, we are solving the complex kinetic 
equations simultaneously for two power modulations: 
Py cos (oot) and Py sin (gX) . The solution to Eq. (13) is 
then represented by the coefficient of e^*"'. In the case 
where the response follows the sinusoidal modulation 
without any lag, the response will be given by Nq''"' 
where A'^ is a real number If the response of the system 
lags the modulation then there will be a measurable 
phase lag between the sinusoidal excitation and the 
sinusoidal response. In this case the solution will be 
given by Ne''°' where A'^ is a complex number The 



complex number A'^ can be written as jA^e-''' where |A^ is 
the magnitude and 6 is the phase difference between 
the driving modulation and the complex solutions. The 

phase is calculated using 0= tan"'( — — — ) where Im(7V) 
^ ^ Re(Af) 

is the imaginary part of the complex number and 
Re(A') is the real part. For small values of the phase, 
0=Im(AO/Re(AO(<4% error for 0< 17°). Since the 
phase lag is approximately given by the ratio of the 
imaginary and real parts of the solutions, the calcula- 
tion of the phase is very simple. The mathematical 
economy of this method is the motivation for using a 
complex modulation. Since the modulation contains a 
constant component, P^, the solution is also assumed to 
have a constant component, A'q, and a time dependent 
response, Ny. The explicit assumed solution has a com- 
plex Fourier expansion given by 



{N{t)) = N, + N,e^"' 



(14) 



where we will ignore terms of higher frequencies 
because they do not contribute to the measured fluores- 
cence response. The procedure for finding the solution, 
{N{t)), is to insert Eq. (14) into Eq. (13) and set equal 
the terms on each side that have the same time varia- 
tion. Carrying out the prescribed procedure, the time 
independent part and the part that varies as e^'^are 
given by 



N, 



N,-- 



AN" 



A+rjaPg 
B + j(o 



AN" 

B 
-N^riaP^B 

B'+G)^ 



.N^r]aPo\ 



B'+of 



(15) 



The constant A is defined by ^ = v/1.274w, and the 
constant B is defined hy B=A + rjaP^. According to 
Eq. (12), for very small power levels the integrated 
fluorescence signal is proportional to P{i){N{t)). The 
fluorescence response is obtained by inserting the 
expression for power of the laser, P{t), and the solution 
{N{t)) given by Eq. (15), yielding the expression for the 
complex magnitude of the fluorescence signal: 
N^Py + N^^Pq +jNy^P(,. After some algebra, we obtain 
the ratio of the imaginary and real part of the fluores- 
cence signal given by 



(16) 



imaginary 


_ N,nPo 


rjamPf, 


real 


A^O^+A^Re^O 


A(A+riaPg) + a)^ 
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The form of Eq. (16) suggests a possible interpreta- 
tion of the measurements. An approach is to measure 
the ratio of the imaginary and real parts of the response 
as a function of frequency for several power levels P^. 
At each power level it is possible to fit the frequency 
dependence to a function of the form ayCdliai + of) and 
obtain the value of the parameter Uy. According to 
Eq. (16), the parameter value should depend linearly on 
power with the slope equal to 77a = k^ ha. From the 
slope it is possible to determine the product of the 
decay rate, molecular rate constants, and a property of 
the laser beam. This is not very satisfactory because the 
result does not provide unique information about the 
photodegradation rate. 

Consider next the case of an larger power level. 
Repeating the steps leading to Eq. (15) and using the 
entire Eq. (9) leads to the following solution 



Af„ 



Af, 



■TV* 



A{\ + ha'P^) 
riaPg+Ail+ba'Pg) 

-ria ; — Nf,P. 

i\+baP^)i\ + ba'P^) ° " 



(17) 



riaPf^+A(\+ba'Pf^) 
(l+baP.) 



+ J(0 



Applying this solution to calculate the fluorescence 
response as given by Eq. (12), we obtain the ratio of the 
imaginary part and the real part given by 



riaP^ 



imaginary 



(\+baP^)(\+ba'P^) 



-03 



real rjaP^+Ajl + ba'P^) 
(l + baPg) 



A{\ + ba'P,) + o)^ 



(18) 



At first glance this result does not look promising. 
However, the denominator is to a good approximation 
dominated by the term A'^ + of and the ratio of a and a' 
is approximately a constant for a given laser beam 
geometry. This ratio can be calculated by Eq. (10). 
Figure 3 shows the value of a and c/ as a function of 
beam width for three values of the photodegradation 
rate, k^. The beam width is one half of the beam waist. 
In the first approximation, the quantities a and c/ are 
related to each other by a constant factor of about 1.45. 
For specific beam width of 20 nm the value of the ratio 
ranges from 1.42 to 1.66 depending on the value of the 
photodegradation rate constant. In absolute terms, a 
range of 300 % variation in the photodegradation rate 
constant results in a 14 % variation in the ratio a'/a. 




lb 2D 

beam width, |am 



Fig. 3. The dependence of the ratio a'/a on the beam width and the photodegradation rate constant. The three lines (dashed, dotted, and soUd) 
correspond to the indicated values of the photodegradation rate. The calculation is performed using the steady state solutions of the kinetic model. 
The calculation shows that the ratio, a'/a, varies relatively little when photodegradation rate varies by a factor of four 
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Thus it is reasonable to set the value of the ratio aVa 
to a constant in Eq. (18) and proceed to evaluate the 
photodegradation rate constant. This is demonstrated in 
the following. 

3b. Proposed Measurement Procedure 

Equation (18) suggests that the ratio of quadrature to 
the in-phase components of the measured response can 
be fitted to a function of modulation frequency of the 
form 



imaginary 



a,(0 



real a^+co 



(19) 



leaving only the photodegradation rate. The final result 
is summarized by 



*-f 



(21) 



A practical consequence of Eqs. (20) and (21) is that 
the absolute power level of the laser beam does not 
have to be measured. Using the absolute power level or 
the relative power in Eq. (20) simply changes the mag- 
nitude of the two parameters c and d, however the ratio 
of the two parameters remains the same in both cases. 

In the following section, measurements will be 
presented of the photodegradation of erythrosine B in 
order to demonstrate the validity of the proposed 
measurement technique. 



Measurements carried out at different values of 
power provide a dependence of the parameter ai on 
laser power. Equation (18) suggests that the measured 
parameter a^ obtained from the frequency dependence 
of the response should depend quadratically on the 
laser power. Hence it should be possible to extract two 
parameters from the fit of the dependence of ai on 
power Setting the ratio a' I a equal to 1 .44, the explicit 
functional form of the dependence of ai on power is 
given by 



cP, 



{\ + dP^){\+d\AAP^) 



(20) 



The two parameters are c = rja as in the case of 
low power, and the additional parameter d=ba. The 
ratio of these two parameters gives the photodegra- 
dation rate. What is the origin of this happy fact? The 
variation of a i on laser power has two sources, first due 
to the photodegradation process described by para- 
meter c, and second due to the saturation of the excited 
triplet state population described by parameter d. 
The reduction due to photodegradation involves the 
same molecular optical parameters as the reduction due 
to saturation. Therefore in the ratio of the two fit 
parameters, eld, the optical parameters cancel out 



4. Measurement of the Photodegradation 
Rate of Erythrosine B 

Figure 4a shows the chemical structure of erythrosine 
B, a derivative of fluorescein with four hydrogen atoms 
replaced by iodine atoms. The heavier iodine atoms 
increase the rate of transitions from the excited singlet 
state to the slightly lower first excited triplet state. 
Recent measurements have shown that erythrosine B 
exhibits phosphorescence at 700 nm (relaxation of the 
triplet state), and delayed fluorescence which occurs 
due to backward transitions from the triplet to the sin- 
glet excited state and subsequent radiative relaxation of 
the singlet state [6]. The emission spectrum, shown in 
Fig 4b, did not show significant intensity in the 700 nm 
region. Measurements were performed below 25 °C to 
minimize delayed fluorescence which becomes signifi- 
cant only at elevated temperature (> 30 °C) [7]. The flu- 
orescence intensity from erythrosine B depends non 
linearly on the intensity of the incident illumination. 
This dependence was interpreted as a depletion of the 
ground state population due to the saturation of the 
triplet excited state [7-9]. The large triplet state popula- 
tion and the quenching of the triplet state by dissolved 
oxygen has been utilized in the design of dissolved 
oxygen sensors [10]. Because of the large triplet state 
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(a) 



(b) 




wavelength, nm 

Fig. 4. (a) The chemical structure of erythrosine B. The fluorophore 
is obtained by replacing four hydrogen atoms in fluorescein by 
iodine atoms. The heavy iodine atoms increase the transition rate 
between the singlet and triplet excited states, (b) The normalized 
absorption and emission spectra of erythrosine B in phosphate buffer 
saline (PBS). The fluorescence emission is insignificant beyond the 
wavelength of 620 nm. 

population, erythrosine B is a good candidate for demon- 
strating the fluorescence modulation technique. 

A schematic of the measurement apparatus is shown 
in Fig. 5 which shows the apparatus from the top. A 
diode laser (Ventus 532 System, Laser Quantum Ltd, 
Stockport, UK)' serves as the source of illumination. The 
laser beam passes through a filter wheel containing a set 
of neutral density (ND) filters with the optical density 
(OD) varying from to 1 in increments of 0.2. Turning 
the filter wheel provides a means to reproducibly vary 
the laser power The laser itself was run at a indicated 



Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 



current of 60 mA corresponding to a nominal power 
output of about 30 mW. After passing the filter wheel, 
the laser beam strikes a mechanical shutter 
(CH-60, Boston Electronics, Brookline, MA) that 
provides the modulation of the laser beam intensity. 
Finally, the laser beam is focused by a 4X objective in a 
flowing solution contained in a flow thru cuvette. The 
cuvette flow channel has a cross section of 1 mm by 
4 mm, and the laser beam propagates in the z direction 
defined as that parallel to the 10 mm side of the cuvette. 
The flow is out of the plane of Fig. 5 (defined as the x 
direction). The photomultiplier detector (PMT), placed 
along thej' direction, views the illuminated portion of the 
solution through a lens, a band pass filter, and two aper- 
tures. The first aperture controls the amount of light 
collected and the second aperture (placed at the focal 
plane) determines the part of the illuminated solution 
which is observed. The band pass filter, centered at 
560 nm, ensures that only erythrosine B fluorescence is 
detected. The signal from the PMT is routed to a current 
amplifier (SR570, Stanford Research Systems, Stanford, 
CA) and finally to a lock in amplifier (SR 830 DSP 
Lock-In, Stanford Research Systems, Stanford, CA). The 
lock-in amplifier is controlled by a computer which sets 
the lock-in internal reference frequency and reads the 
resulting PMT signal. A TTL output, locked to the inter- 
nal reference signal of the lock- in amplifier and located 
on the rear panel of the lock-in amplifier, is used to ener- 
gize the shutter controller The time constant of the lock- 
in amplifier was set to 1 s. During the acquisition of the 
signal, a time delay was placed between the setting of the 
lock- in reference frequency and the reading of the PMT 
signal. The delay was set to the time it took for the lock- 
in amplifier to settle after a change of reference frequen- 
cy. The acquisition of signal consisted of 20 sequential 
"reads" of the in-phase and quadrature components. The 
ratio of the quadrature and in-phase components was 
calculated and presented as the phase. Thus the phase 
presented in the following figures is the average of the 20 
"reads". The standard deviation of the 20 samples was 
much smaller than the size of the circles in the plots 
presented in Figs. 6a and 7a. Twelve frequencies were 
chosen such that the logarithm of the frequencies was 
spaced equally. This resulted in more points at low 
frequencies where the response variation was the great- 
est. A 35 mL capacity syringe pump (Harvard Apparatus 
22, South Natick, MA) was used to pump the solution 
through the cuvette. At the low flow rates used for the 
measurements, six frequency scans were performed each 
for a different ND filter wheel setting. A constant solu- 
tion flow was maintained throughout the six frequency 
scans. 
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Fig. 5. A schematic of the apparatus used for the measurement of the photodegradation rate. The flow thru cuvette is viewed from the top with 
the flow coming out of the paper (x axis in Fig. 1). The laser beam passes thru the cuvette along the long dimension of the flow channel. The 
filter wheel contains ND filters of increasing CD, and the shutter modulated the laser intensity with a ftequency controlled hy the lock-in ampli- 
fier internal reference oscillator. The signal detected by the PMT is routed to the lock-in amplifier which gives the in-phase and quadrature 
components of the signal. 



The open circles in Fig. 6a show a typical response for 
erythrosine B solution for a filter wheel set to OD = 0. 
The solid circles in Fig. 6a show the response from a 
sintered plate of glass replacing the fluorophore solu- 
tion. The scattering from glass is an instantaneous 
process and it should not lead to any phase shift in the 
observed signal relative to the excitation. Therefore the 
observed phase shift is due to the behavior of the instru- 
mentation. The shutter has a fixed opening and closing 
time interval, therefore the shape of the light pulse will 



change with frequency resulting in an instrumental 
"apparenf phase shift. The instrumental phase shift 
was highly reproducible, and it was always subtracted 
from the phase shift observed for the fluorophore 
solution, resulting in a phase shift which depended only 
on the fluorophore population dynamics. As an exam- 
ple, the phase shift resulting from the subtraction of the 
two measured phase shifts in Fig. 6a is shown in Fig. 6b 
by the solid circles. The solid curve in Fig. 6b is a 
spline curve used to show the trend in the data. 
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Fig. 6. (a) The solid circles show the ratio of the quadrature component and the in-phase component of the 
measured scattered signal as a function of the modulation frequency. Since scattering is an instantaneous event 
the measured phase is due to the response of the instrument. The open circles show the ratio of the quadrature 
component and the in-phase component of the measured fluorescence signal as a function of the modulation 
frequency. The fluorescence response is delayed relative to the illumination modulation resulting in additional 
phase shift, (b) Subtracting the instrumental phase shift (soUd circles in Fig. 6a) from the fluorescence phase 
shift (open circles Fig. 6a) gives the phase shift due to the dynamics of the fluorescence process alone. The solid 
curve in Fig. 6b is a spHne fit to the data showing the trends in the data. 
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Figure 7a shows a typical set of six frequency scans 
taken at decreasing laser power for erythrosine B in a 
PBS solution. The instrumental phase shift was sub- 
tracted out in all cases. The responses in Fig. 7a were 
fitted to Eq. (9) and the value of parameter Uy was 
obtained as a function of laser power (see Eq. (20). The 
solid circles in Fig. 7b show the dependence of the 



parameter a^ obtained from the data in Fig. 7a on laser 
power. The uncertainties in the values of the parameters 
are of the size of the solid circles and are indicated by 
the error bars. The uncertainty estimates were obtained 
from the fitting algorithm used in the SigmaPlot 9 
application (Systat Software, Inc., San Jose CA). The 
solid curve in Fig. 7b is a fit to Eq. (20). The deviations 
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Fig. 7. (a) The six sets of symbols give the response at six different laser power levels. The response at each power level is fitted to the expect- 
ed response given by a^iaKa^ + la'') to obtain the parameter a,, (b) The solid circles show the dependence of a, on the relative power level of the 
incident laser light. The error bars associated with the values of Oj were obtained from the fitting procedure used in 
Fig. 7a. The solid curve is a best fit to Eq. (20) discussed in the text. The ratio of the fit parameters gives the photodegradation rate. 
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between the data and the fit are small (< 5 %) yielding 
good estimates of the parameters c and d in Eq. (20). 
The ratio of the two parameters gives a photodegrada- 
tion constant of (33 ± 2) s~' where the uncertainty is 
discussed in Sec. 5. An important point is that the 
power axis in Fig. 7b is given in terms of relative power 
obtained from the value of the OD filter Using the 
absolute power levels would change the values of the 
fit parameters but the ratio of the parameters would 
remain the same. This suggests that the value of the 
absolute power of the laser beam is not critical in 
obtaining the photodegradation rate. There are system- 
atic errors in the analysis inherent in the assumption 
that the ratio a' I a is a constant. The present analysis 
suggests that the systematic errors are relatively small. 
The measurements were repeated in solutions of 
increasing viscosity. The viscosity was progressively 
increased by mixing in 1 % by mass dextran (a poly- 
saccharide) or 1 % by mass PEO (polyethylene oxide) 
long chain polymers in the buffer containing erythro- 
sine B. The response and analysis of the fluorescence 
measurements were similar to those shown in Fig. 7a 
and Fig. 7b. Table 1 shows the measured photodegrada- 
tion rates in the increasingly viscous solutions. The fact 
that the photodegradation rate decreases with increas- 
ing viscosity suggests that the underlying reaction is 
diffusion limited. This would be the case if erythrosine 
B reacted with the photoproduct resulting from the 
relaxation of the triplet state. For example singlet 
oxygen can be produced by a reaction between the 
dissolved triplet oxygen and the excited triplet state of 
erythrosine B. Once produced, the singlet oxygen can 
react with erythrosine B and alter its structure. The role 
of diffusion in the chemical reaction is on the micro- 
scopic level and does not affect the assumptions made 
previously that the role of diffusion is negligible in 
macroscopic mass transfer into the entire illuminated 
region. The dependence of photodegradation rate on 
viscosity (and hence on diffusion) suggests the 
approach can be used to study the mechanism of the 
photodegradation process. The photodegradation rate 
can be measured for a set of environmental conditions 
that are deemed important. The variation of the pho- 
todegradation rate with the environmental parameters 
can suggest a model for the photodegradation process. 
This strategy was followed previously by relating the 
decrease in the photodegradation rate of fluorescein to 
the increase of the concentration of «-propyl gallate, a 
known singlet oxygen quencher [11]. 



Table 1. 



Solvent 

PBS buffer 

PBS buffer + 1 % dextran 

PBS buffer + 1 % PEO 



kj=c Id, s 
33 ±2 
28 ±2 
20 ±2 



5. Uncertainty Analysis 

The photodegradation rate is given by the ratio of 
two parameters, c and d, in Eq. (20). The parameters 
are obtained by fitting the function containing the two 
parameters to the values of the parameter ai obtained 
for different values of the relative laser power The 
values of ai are in turn obtained by fitting the function 
containing a^ Eq. (19) to the measured ratios of quadra- 
ture and in phase fluorescence signals at different 
modulation frequencies. The uncertainties of the two 
measured fluorescence signals are set equal to the 
standard deviation of several repetitive measurements 
taken under identical conditions. Given the uncertain- 
ties of the quadrature and in phase signals, the uncer- 
tainty in the ratio of the two signals can be obtained by 
using a simple error propagation formula. In theory, the 
best fit to the ratio of the signals, using the uncertain- 
ties as weights, gives values of ai with uncertainties 
which take into account both the uncertainties in the 
ratio of the signals and the fitting process. The uncer- 
tainties in Qi would in turn be used as weights in the 
fitting to obtain the parameters c and d. In theory, the 
final uncertainties of c and d would depend on the 
uncertainties of the quadrature and in phase signals plus 
the additional uncertainties introduced by the two 
fitting procedures. In practice, the SigmaPlot algorithm 
used in the fitting does not use uncertainties of the data 
in the fitting algorithm. Rather the algorithm estimates 
the parameter uncertainties by examining the devia- 
tions between the unweighted data and the best func- 
tional fit. The uncertainties presented in Table 1 are the 
uncertainties given by the fitting algorithm used to 
obtain the parameters c and d. 



6. Conclusion 

The fluorescence response of fluorophores in solu- 
tion illuminated by a modulated laser beam can be 
measured using synchronous detection with a lock-in 
amplifier The ratio of the quadrature to the in-phase 
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components of the measured response can be fitted 
to a function of the modulation fi-equency given by 
ai(ol{a2+ of). The fitting yields values of the para- 
meters fli and flj- In this paper we show how the analyt- 
ical solution of a mathematical model of the experiment 
(derived in an earlier paper [3]) allows the expression 
of the parameters Ui and aj in terms of the physical 
constants of the experiment. The measurement can be 
repeated for several power levels resulting in a meas- 
ured dependence of the parameter Oi on the power of 
the laser. The dependence of the parameter Oi on the 
power can be fitted to a function, ai = cPQ/{\ + dP^ 
(1 +d\.44P^, where 1.44 is a property of the laser 
beam geometry corresponding to the beam width. The 
ratio of the two fitting parameters, c Id, gives the pho- 
todegradation rate k^. The ratio is not sensitive to 
uncertainties in the absolute power level of the laser. 
The utility of this measurement method is based on the 
separation of the photodegradation rate from all of the 
optical transitions rates associated with the molecular 
excited states. 
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